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IRATORY  POSITION  STATEMENT 


This  report  documents  an  engineering  analysis  and  the  resulting  computer 
program  which  are  considered  to  be  technically  sound  for  the  pre- 
diction of  helicopter  fuselage  flow  fields  without  rotor  wake  effects. 
The  computed  results  from  the  program  should  be  used  for  design  purposes 
only  after  correlation  of  the  program  with  test  data  for  an  aircraft  of 
similar  configuration. 


The  computer  program  resulting  from  this  contract  will  be  provided  to 
qualified  users,  upon  request,  for  use  In  the  design  and  analysis  of 
rotary-wing  aircraft. 


The  project  engineer  for  this  contract  was  Mr.  W.  D.  Vann,  Aeromechanics 
Technical  Area,  Aeronautical  Technology  Division. 
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. aSITHACT  (Continue  on  nnrH  ■(<*•  II  nococoocy  ond  Identify  by  block  mmboc) 

N A computer  program  has  been  developed  which  allows  the  flow  around 
general  aircraft  configurations  to  be  modeled,  including  the  effects  of 
regions  of  separated  flow.  The  program,  although  developed  specifically  to 
handle  the  regions  of  separated  flow  associated  with  the  typically  bluff 
helicopter  airframes,  retains  the  capability  to  model  the  more  slender  shapes 
associated  with  fixed-wing  aircraft.  The  flow  around  the  specified  configu- 
ration is  calculated  iteratively,  the  first  step  being  the  solution  of  the  — 
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basic  potential  flow.  This  is  followed  by  the  definition  of  the  potential 
flow  streamlines  over  the  surface,  the  calculation  of  boundary  layer  growth 
along  the  streamlines,  the  identification  of  the  separated  flow  region  and 
the  eventual  recalculation  of  the  potential  flow,  including  the  effects  of 
boundary  layer  growth  and  separated  flow.  The  iterations  between  potential 
and  viscous  flow  may  be  repeated  as  often  as  desired.  During  its  development, 
the  program  has  been  exercised  on  a number  of  helicopter  and  aircraft  con- 
figurations. This  experience  has  shown  that  the  program  can  be  used  to  model 
the  flow  about  the  involved  shapes  typical  of  the  new  generation  of  rotor- 
craft  now  going  into  service  without  any  increase  in  difficulty,  either  in 
preparing  the  input  data  or  in  executing  the  solution,  when  compared  with 
earlier  purely  potential  flow  solutions,  tv 
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This  program  was  sponsored  by  the  Applied  Technology  Laboratory,  U.S.  Army 
Research  and  Technology  Laboratories  (AVRADCOM) , Fort  Eustis,  Virginia,  and 
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INTRODUCTION 


Until  very  recently,  design  for  low  drag  has  never  been  a major  constraint  in 
the  development  of  rotary  wing  aircraft  since  with  cruise  speeds  never  much 
more  them  100  knots,  the  contribution  of  the  airframe  to  the  total  power  con- 
sumed has  been  relatively  minor  when  compared  with  the  rotor  effects.  However, 
with  the  cruise  requirement  raised  to  145  knots,  as  it  has  been  for  both  the 
Armed  Attack  Helicopter  (AAH)  and  the  Utility  Tactical  Transport  Aircraft 
System  (UTTAS),  airframe  drag  is  playing  a more  significant  part.  This,  to- 
gether with  configuration  design  constraints  introduced  by  other  factors  (such 
as  air  transportability ) and  the  need  to  define  achievable  performance  relative 
to  a very  challenging  specification,  underscored  the  need  for  a reliable  means 
of  predicting  drag  during  the  design  process. 

Traditional  empirical  drag  estimation  methods  are  notoriously  imprecise,  even 
when  used  within  a given  generation  of  rotorcraft,  and  when  used  to  evaluate 
a design  where  there  is  a jump  in  the  level  of  technology,  they  cannot  be 
expected  to  be  reliable.  Where  the  configuration  is  slender,  such  as  most 
fixed  wing  aircraft  in  the  cruise  configuration,  and  where  there  are  no  regions 
of  separated  flow,  potential  flow  solutions  and  a subsequent  boundary  layer 
and  skin  friction  calculation  provide  a good  estimate  of  drag.  Unfortunately, 
on  the  typical  helicopter  configuration  bluff  shapes  are  the  rule  rather  them 
the  exception,  regions  of  separated  flow  are  present  on  most  of  the  aft- facing 
surfaces,  and  potential  flow  methods  are  inadequate.  It  was  to  provide  a 
method  by  which  paneling  methods  developed  for  potential  flow  field  modeling 
could  be  extended  to  the  bluff  body  helicopter  problem  that  the  work  described 
in  Reference  1 was  carried  out. 

In  the  study  described  in  Reference  1,  a potential  flow  method  in  which  the 
body  was  represented  by  a number  of  flat  source  panels  was  used  as  the  basis 
for  a streamline  and  boundary  layer  calculation.  Separation  location  was 
identified  on  each  streamline,  and  with  the  presence  of  the  boundary  layer 
thickness  and  separated  flow  region  modeled,  the  calculation  was  recycled 
through  the  flow  field  calculation.  This  process  was  repeated  as  often  as 
was  needed  to  converge  on  a solution. 

The  computer  program,  code  named  DRAG,  demonstrated  the  ability  to  predict  as 
part  of  the  aerodynamic  calculation  the  base  pressure  in  separated  flow 
regions.  The  program  was  deficient  in  two  areas,  however.  First,  in  order 
to  take  advantage  of  an  existing  streamline  tracing  procedure,  each  column 
of  panels  was  constrained  to  have  the  same  number  of  rows.  This  assumption 
of  a regular  panel  grid  resulted  in  considerable  reduction  of  streamline 
computation  effort,  but  made  the  task  of  paneling  the  configuration  somewhat 
more  difficult,  the  user  being  restricted  to  simple,  regular  bodies,  or  if 


1.  Dvorak,  F.A.,  Maskew,  B.,  and  Woodward,  F.A.,  "Investigation  of  Three- 

Dimensional  Flow  Separation  on  Fuselage  Configurations",  USAAMRDL-TR-77-4, 
Eustis  Directorate,  U.S.  Army  Air  Mobility  R & D Laboratory,  Fort  Eustis, 
Virginia,  March  1977,  AD  A039382. 
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complex  bodies  were  required,  to  inefficient  panel  densities  in  some  regions 
of  the  body. 


Second,  close  to  the  separation  point,  there  existed  a sharp  excursion  in 
the  calculated  pressures  which  was  absent  in  the  experimental  data.  The 
reason  for  the  anomaly  lies  in  the  type  of  singularity  chosen  to  model  the 
body.  In  this  case,  carrying  over  from  an  earlier  model,  constant  source 
panels  were  used  on  the  body  surface,  while  the  separated  wake  was  represented 
by  constant  vorticity  panels.  Consequently,  an  abrupt  change  in  both  singu- 
larity type  and  strength  was  introduced. 

The  present  study  was  undertaken  to  remove  the  constraints  existing  in  the 
DRAG  program.  As  part  of  this  effort,  several  improvements  and  additional 
capabilities  were  provided  for  the  user  of  the  program.  The  patching/paneling 
restriction  was  removed  by  going  to  a multi-patch  representation  of  the  body, 
retaining  within  each  patch  the  regular  structure  required  by  program  DRAG 
while  requiring  only  the  patch  joint  information  in  addition  to  the  basic  input. 
This  provides  a very  efficient  and  flexible  paneling  scheme.  The  irregularity 
in  the  pressures  calculated  was  removed  by  going  to  a higher-order  singularity 
in  place  of  the  simple  constant  strength  sources  used  in  the  original  code. 
Additional  features  added  include  the  ability  to  model  wind  tunnel  walls  and 
floor,  and  engine  efflux  modeling. 


The  program  in  its  present  form  has  been  used  to  model  the  flow  about  represen- 
tative helicopter  configurations/  and  examples  of  its  use  are  presented  in  the 
report . 
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GEOMETRY  DEFINITION 

As  in  the  original  version  of  the  DRAG  program,  a configuration  is  rep- 
resented geometrically  by  a network  of  interconnected  panels.  Within 
the  program  there  are  a number  of  routines  which  require  information 
from  neighboring  panels  (e.g.,  streamline  tracing,  boundary  layer  source 
and  skin  friction  redistribution,  surface  vorticity  continuity) . This 
information  is  readily  available  from  a rectangular  grid  of  panels  as 
used  in  the  basic  version  of  the  DRAG  program.  A single  rectangular 
grid,  however,  is  very  restrictive  for  applications  to  practical  prob- 
lems and  so  a scheme  has  been  developed  based  on  multiple  PATCHES  of 
peine  Is  (Figure  1) . Each  patch  has  a rectangular  grid  of  panels  and  so 
maintains  the  ease  of  obtaining  neighboring  panel  information  within 
itself.  A special  SEARCH  routine  (described  in  a later  section)  has 
been  developed  to  determine  information  on  neighboring  panels  across 
the  patch  edges  where  there  is  not  necessarily  a direct  correlation 
between  neighboring  panels.  Each  patch,  therefore,  is  a complete  unit 
in  the  sense  that  its  panel  spacing  is  not  necessarily  related  to  a 
neighboring  patch.  In  this  way,  very  general  configurations  can  be 
paneled.  Clearly,  however,  the  number  of  patches  should  be  kept  as 
small  as  possible,  consistent  with  reasonable  panel  shapes.  Patch 
surfaces  (when  opened  out)  should  be  usually  four-sided,  although  one 
side  or  two  opposite  sides  can  degenerate  to  points,  producing  tri- 
angular panels  locally. 

The  patches  are  numbered  automatically  by  the  program  as  it  constructs 
the  peine 1 model  in  WBPANX,  the  patch  number  being  a function  of  the 
order  in  which  the  section  input  is  read.  The  key  parameter  in  deter- 
mining patching  is  the  parameter.  Flag,  on  Card  10B  of  the  input  deck. 
When  a value  of  1 is  read  for  this  parameter,  the  code  knows  that  the 
last  section  of  a patch  is  being  read.  It  then  assigns  a patch  number 
to  the  panels  constructed  since  the  last  section,  read  with  FLAG  = 1, 
records  the  number  of  rows  and  columns  in  the  patch  and  the  number  of 
the  last  panel  in  the  patch.  This  information  is  used  in  the  solution 
of  the  influence  coefficient  matrix  in  the  streamline  calculation. 

The  number  (order)  of  the  patches  may  be  deduced  from  the  input  data 
deck  and  has  to  be  known  in  order  to  supply  the  streamline  code  with 
information  on  neighboring  panels.  This  is  supplied  in  the  form  of 
a table  of  patch  side  overlaps,  the  patch  sides  being  numbered  anti- 
clockwise when  viewed  from  outside  the  body  with  side  one  being  that 
side  closest  to  the  front  of  the  body.  A typical  patch  schematic 
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for  the  body  shown  in  Figure  1 is  given  in  Figure  2.  A partial  list  of 
the  patch  joint  information  is  given  below. 
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Since  the  body  is  symmetric,  the  joint  along  the  middle  of  the  model  is 
a reflection  with  the  patch  being  overlapped  by  itself  (on  the  other 
side).  Hence  1,4  overlaps  1,4.  Also,  even  though  side  1 of  patch  3 is 
a point  (actually  two  points  superimposed) , it  must  be  modeled  if  pas- 
sage of  a streamline  through  that  region  is  desired. 


POTENTIAL  FLOW  MODEL 


A new  potential  flow  model  has  been  developed  to  replace  the  original  con- 
stant source  panel  representation  described  in  Reference  2.  The  new  aero- 
dynamic singularity  consists  of  a combination  of  constant  and  linearly 
varying  vortex  distributions  on  the  body  panels.  This  type  of  singularity 
has  been  successfully  applied  to  the  analysis  of  separated  flows  on  two- 
dimensional  airfoils  (Ref.  3).  A generalization  of  these  vortex  singu- 
larities suitable  for  the  analysis  of  arbitrary  body  shapes  is  derived 
below . 

Consider  a typical  body  panel  having  unswept  leading  and  trailing  edges, 
arbitrary  edge  slopes,  m^  and  m2,  and  chord  length,  c. 


For  convenience,  the  origin  of  coordinates  is  located  at  the  intersection 
of  the  leading  edge  and  the  inboard  side  edge.  The  aerodynamic  influence 
of  the  panel  can  be  obtained  by  integrating  the  influence  of  an  elementary 
vortex  filament  and  its  wake,  from  the  leading  edge  to  the  trailing  edge. 

Considering  the  contribution  of  the  elementary  bound  vortex,  the  induced 
perturbation  velocity  components  may  be  written 
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2.  Woodward,  F.A.,  Dvorak,  F. A. , and  Geller,  E.W.,  "A  Computer  Program 
for  Three-Dimensional  Lifting  Bodies  in  Subsonic  Inviscid  Flow", 
USAAMRDL  TR-74-18,  Eustis  Directorate,  U.S.  Army  Air  Mobility  Research 
and  Development  Laboratory,  Fort  Eustis,  Va.,  April  1974,  AD  782202. 

3.  Maskew,  B.,  and  Dvorak,  F.A.,  "The  Prediction  of  Cg,  Using  a Sepa- 
rated Flow  Model",  J.  Am.  Hel.  Soc. , April  1978.  ^"aX 
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r2  = (mx  - y)2  + (1  + m2)z2 


and  k = Q for  linear  edge  vortex  associated  with  constant  vortex  dis- 

tribution i and 

k = 1 for  quadratic  edge  vortex  associated  with  linear  vortex 
distribution. 


Evaluating  the  above  integrals , the  following  results  are  obtained  for 
the  three  components  of  velocity  at  a point  (x,y,z)  induced  by  the 
bound  vortex  distributions. 


Constant  Bound  Vortex  Distribution 


Linear  Bound  Vortex  Distribution 


These  expressions  give  the  influence  of  a single  corner  of  the  panel,  with 
the  origin  of  coordinates  located  at  that  corner.  The  influence  of  the 
complete  panel  is  obtained  by  summing  the  influences  of  each  of  the  four 
corners . 


AERODYNAMIC  REPRESENTATION 

The  body  surface  is  subdivided  into  a number  of  patches,  each  patch 
having  a fixed  number  of  rows  and  columns.  Within  each  patch,  the  con- 
stant and  linearly  varying  vortex  singularities  are  formed  into  over- 
lapping wedges  over  pairs  of  adjacent  panels  along  each  column,  as  il- 
lustrated below. 


Vortex  Wedge  Trailing 
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The  aerodynamic  influence  of  each  vortex  wedge  includes  the  contributions 
of  the  constant  and  linearly  varying  bound  vorticity,  the  edge  vorticity/ 
and  a pair  of  concentrated  vortices  trailing  downstream  in  the  wake.  The 
unknown  vortex  strength  is  the  height  of  the  wedge,  while  the  control 
points  (used  for  satisfying  the  boundary  conditions  of  tangential  flow) 
are  located  at  panel  centroids. 

For  convenience  in  accounting  for  the  influence  of  the  vortex  wake, 
adjacent  pairs  of  vortex  wedges  are  combined  to  form  doublets.  The 
strengths  of  the  vortex  wedges  are  weighted  according  to  the  chord 
lengths  of  the  panels  such  that  the  trailing  vortices  from  each  wedge 
exactly  cancel  in  the  wake.  The  resulting  doublet  singularity  extends 
over  three  adjacent  panels  as  shown  below.  The  control  point  is  loca- 
ted at  the  centroid  of  the  middle  panel. 


Combined  Doublet 
Singularity 
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In  the  application  of  the  doublet  singularity,  special  care  is  required 
at  the  junction  of  two  patches  to  ensure  conservation  of  vorticity.  The 
last  two  rows  on  the  upstream  patch  contain  truncated  doublet  singulari- 
ties. Vorticity  is  conserved  within  this  patch  by  introducing  stepwise 
constant  line  vortices  along  the  trailing  edge  of  the  patch.  The 
strengths  of  these  edge  vortices  are  proportional  to  that  of  the  corre- 
sponding truncated  doublet.  This  concentrated  trailing-edge  vorticity  is 
then  transferred  to  the  leading  edge  of  the  downstream  patch.  This  in 
turn  is  fed  into  the  nearest  side  edges  of  the  first  row  of  panels  in 
this  patch  and  distributed  linearly  over  each  of  these  panels.  In  this 
way  any  residual  vorticity  can  be  exactly  cancelled  by  the  first  row  of 
panels.  The  influence  of  the  remaining  doublet  singularities  in  the 
downstream  patch  can  then  be  computed  in  the  usual  way. 


POTENTIAL  PLOW  SOLUTION 

The  doublet  singularities  are  used  to  calculate  the  normal  velocity  at 
panel  control  points,  resulting  in  the  formation  of  a large  system  of 
linear  equations  relating  the  doublet  strengths  to  the  prescribed  normal 
velocity  of  each  control  point.  However,  the  use  of  doublet  singulari- 
ties alone  is  not  sufficient  to  ensure  a unique  solution.  A constant 
source  singularity  is  also  distributed  over  each  panel  and  used  to  cal- 
culate the  normal  velocity  at  the  control  points  around  the  first  ring 
of  panels  at  the  body  nose.  The  resulting  system  of  equations  may  then 
be  solved  to  determine  the  unknown  source  and  doublet  strengths. 

The  ntaerical  stability  of  this  system  of  equations  can  be  improved  by 
introducing  additional  constant  source  singularities  on  each  panel.  The 
strengths  of  these  additional  sources  are  set  equal  to  the  axial  com- 
ponent of  the  onset  flow.  For  onset  flows  parallel  to  the  x-axis,  this 
is  equivalent  to  applying  the  Green's  Identity  formulation  described  in 
Reference  4.  This  additional  source  distribution  has  the  effect  of  re- 
ducing the  doublet  strengths  required  to  satisfy  the  boundary  condition 
equations,  and  significantly  improves  the  numerical  stability  and  con- 
vergence of  the  solution: 

The  resulting  system  of  equations  may  be  written  in  matrix  form  as 
follows: 


B.  . a.  + n. 
ij  3 i 


where  A 


ij 


n • u + n • v +n 

"i  Du  y±  Dij 


" w 

Zx  Dij 


B.  . * n • u +n  • V + n • w 

13  -i  Vt  ■>„  \ \j 
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n /4tt  * prescribed  source  strength; 

Xj 


y = unknown  doublet  strength; 


and 


n , n , n - direction  cosines  of  normal  to  control 
xi  zi  point  of  panel,  i; 


u f v , w * three  components  of  velocity  induced  by 
Dij  Dij  Dij  source,  j,  on  control  point,  i»  and 


o,  8 = angles  of  attack  and  yaw  of  onset  flow. 


Once  the  doublet  strengths  have  been  determined,  the  three  components  of 
velocity  at  control  point,  i,  can  be  computed  as  follows: 
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STREAMLINE  CALCULATION 

The  DRAG  program  computes  an  approximate  solution  to  the  problem  of 
three-dimensional  potential  flow  over  an  arbitrary  body.  As  described 
earlier,  the  body  surface  is  now  represented  by  a set  of  patches. 

Within  each  patch  a set  of  plane  quadrilaterals  are  arranged  in  a rec- 
tangular grid;  however,  the  grid  of  panels  can  vary  from  one  patch  to 
another.  The  three  components  of  the  velocity  vector  are  computed  at 
the  panel  centroids.  The  streamline  program  uses  the  velocity  data  at 
four  neighboring  panels  to  provide  a linear  approximation  to  the  veloc- 
ity anywhere  on  a given  panel.  A straight-line  approximation  to  the 
actual  streamline  passing  through  a point  on  the  panel  is  constructed 
and  continued  over  all  panels  upstream  and  downstream.  This  results  in 
relatively  smooth  streamlines  over  the  body,  and  allows  the  geodesic 
curvatures  of  the  streamlines  and  equipotential  lines  to  be  computed. 
From  this  information,  the  metric  or  effective  radius  of  the  body  along 
the  streamline  is  obtained.  The  method  used  follows  closely  that  de- 
scribed in  Reference  5. 


COMPUTATION  OF  THE  PARTIAL  DERIVATIVES 


e 

on  the  f 
fied  by 
stream  b 


XR,YR 


5.  Dawson,  C.W. , and  Dean,  J.S.,  "The  Calculation  of  Streamline  Data  for 
Boundary  Layer  Input",  NSRDC  CMD-4-74,  January  1974. 


If  panel  Q borders  on  a patch  edge,  then  SEARCH  is  called,  and  an  array 
identifying  all  panels  adjacent  to  the  panel  edge  is  generated.  In  the 
example  sketched  below,  the  "shaded"  panels  would  be  listed  as  down- 
stream neighbors  of  panel  Q. 


PATCH  1 


PATCH  2 


PATCH  3 


Initially,  any  one  of  these  panels  is  selected  as  the  downstream  neigh- 
bor, R,  depending  on  the  order  of  description  of  the  joints  between  the 
participating  patches.  Ultimately,  neighbor  R is  selected  by  requiring 
that  the  predicted  intersection  point  be  adjacent  to  the  panel  edge. 
This  procedure  is  explained  in  more  detail  in  the  section  describing 
SEARCH.  The  new  streamline  code  also  ensures  that  the  existing  panel 
(previous  panel  Q)  is  employed  as  a neighbor  to  panel  Q.  Once  the 
neighbors  are  identified,  the  neighboring  panels  are  assumed  to  be  ro- 
tated about  the  common  edges  so  that  all  panels  lie  in  the  same  plane 
as  the  central  quadrilateral.  The  program  then  computes  the  velocity 
and  position  vectors  of  the  neighbors  in  terms  of  the  central  panel  co- 
ordinate system.  From  this  construction,  the  following  partial  deriva- 
tives are  obtained: 
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[30] 


where  VU  is  the  transverse  component  of  velocity  in  the  local  coor- 
dinate system  of  the  upper  neighbor, 

UU  is  the  axial  component  of  velocity  in  the  local  coordinate 
system  of  the  upper  neighbor,  etc., 

DXR  is  the  axial  distance  along  the  surface  between  the  cen- 
troids of  panels  Q and  R, 


DXT  is  the  axial  distance  along  the  surface  between  the  cen- 
troids of  panels  Q and  T, 


DYU  is  the  transverse  distance  along  the  surface  between  the 
centroids  of  panels  Q and  U,  and 

DYL  is  the  transverse  distance  along  the  surface  between  the 
centroids  of  panels  Q and  L. 


The  axial  and  transverse  components  of  velocity  on  panel  Q are  then  ap- 
proximately 
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COMPUTATION  OF  POINTS  ON  A STREAMLINE 

The  streamline  computation  begins  gt  a given  point  in  a specified  quadri- 
lateral with  the  velocity  vector,  V » (U,V) . A stream  function  is  de- 
sired which  will  have  a gonstant  value  of  zero  along  the  streamline. 
Since  the  divergence  of  V is  not  zero,  an  additional  function,  p,  must 
be  found  such  that  the  divergence  of  pv  is  zero.  Then  for  a particle 
following  a streamline, 

d£  _ dy/dt  , V _ £V  (32) 

dx  dx/dt  U pu 
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Thus,  a new  vector  field  is  constructed  whose  streamlines  are  identical 
to  those  of  the  velocity  field,  but  whose  divergence  is  zero.  Hie  func- 
tion, P,  is  given  by  the  series 
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The  resulting  stream  function  is 


SF(xq,yq)  = SFq 
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The  constant  value,  SF0,  is  chosen  so  that  the  stream  function  is  zero 
at  the  specified  point.  Since  the  velocity  is  known  only  through  the 
linear  terms,  it  is  consistent  to  truncate  the  stream  function  after 
the  quadratic  terms. 


Values  of  SF  at  the  four  comer  points  of  the  quadrilateral  are  now  com- 
puted. By  comparing  the  signs  of  SF  at  adjacent  corners,  the  sides 
through  which  the  streamline  passes  are  determined.  If  the  streamline 
passes  through  a side,  the  value  of  SF  is  computed  at  the  midpoint  of 
the  side.  From  these  three  points,  the  intersection  point  is  computed 
from  a three-point  interpolation  formula. 


The  parameter,  t(0  <_  t <_  1),  is  defined  so  that 
x = Xj^  + t(x3  - x^ 

y * yx  + fc(y3  - yx> 
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Here, 
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Then, 

0 = SF(x,y)  = 2SF1(t  - *i)(t  - 1)  - 4SF2t(t  - 1) 
+ 2SF3t(t  - H) 


(36) 


The  root  of  this  equation  is  chosen  so  that  0<  t < 1.  Note  that  if 
there  had  been  two  roots  between  0 and  1,  then  SF  would  have  the  same 
sign  at  both  corner  points,  and  those  intersection  points  would  be  ig- 
nored . 

In  general,  two  intersections  will  be  found  for  the  entire  quadrilateral: 
one  where  the  streamline  enters,  and  one  where  it  leaves.  It  is  possible, 
however,  for  there  to  be  four  intersection  points.  The  next  point  on 
the  streamline  is  chosen  from  the  intersection  points  as  follows:  a 
quantity,  TESTP,  is  computed  for  each  intersection  point  by  taking  the 
dot  product  of  the  vector  from  the  starting  point  to  the  intersection 
point  with  the  velocity  vector  at  the  centroid  of  the  quadrilateral, 
and  then  multiplying  it  by  the  direction  sign.  (The  direction  sign  is 
+1  if  the  streamline  is  being  traced  upstream.)  The  intersection  point 
with  the  maximum  positive  value  of  TESTP  is  chosen.  However/ if  the 
largest  TESTP  is  less  than  or  equal  to  zero,  none  of  the  points  is  ac- 
ceptable. In  general,  this  will  occur  in  places  where  the  panel  spac- 
ing changes  abruptly  across  patch  edges  and/or  where  the  actual  body 
curvature  changes  abruptly.  In  these  cases,  the  approximations  outlined 
above  in  calculating  the  partial  derivatives  can  become  too  large}  con- 
sequently, to  continue  the  streamline,  the  fluid  particle  is  assumed  to 
travel  across  the  panel  from  the  entering  point  at  a velocity  equal  to 
the  velocity  at  the  panel  centroid  as  indicated  in  the  following  sketch. 


Entering  Point 


Panels  bordering  on  a plane  of  symmetry  sometimes  require  special  hand- 
ling. If  the  calculated  intersection  point  lies  on  the  plane  of  sym- 
metry, this  point  is  rejected  and  the  final  intersection  point  is  cal- 
culated geometrically  on  the  upstream  or  downstream  edge  of  finite 
length,  depending  on  the  direction  of  the  streamline  calculation.  The 
sketch  below  illustrates  resultant  streamline  shapes  across  two  dif- 
ferent panels  bordering  on  a plane  of  symmetry. 


If  the  current  panel  borders  on  a patch  edge  and  the  streamline  is  to 
cross  onto  the  next  patch,  then  the  final  intersection  point  is  used  to 
select  the  next  panel  from  the  list  of  adjoining  panels  discussed  earlier. 
Because  there  are  generally  small  gaps  between  the  panels,  this  selection 
process  must  be  carefully  done  to  ensure  continuation  of  the  streamline 
calculation.  For  example,  intersection  points  that  lie  in  an  interstice 
require  special  handling  (see  sketch  below) . 


Possible 
Next  NQ 
Panels 


Here,  an  entering  point  very  close  to  the  corner  point  of  each  possible 
next  panel,  NQ,  is  assumed  in  turn  and  the  streamline  procedure  continues. 
If  the  stream  function  values  on  the  new  panel  indicate  a good  future 
intersection,  then  that  particular  panel  is  chosen  as  the  neighbor  to 
panel  Q in  the  sketch  above.  A final  check  is  made  in  all  cases  to  as- 
sure that  the  final  next  panel,  NQ,  has  actually  participated  in  the 
calculation  of  the  partial  derivatives  associated  with  the  streamline 
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crossing  of  psnel  Q.  If  not,  the  crossing  of  panel  Q is  repeated  with 
NQ  used  as  the  appropriate  neighbor. 

In  sunaary,  the  new  features  provided  in  STREAM  are  the  following: 

1.  SEARCH  identifies  an  array  of  panels  along  a patch  edge 
adjacent  to  panel  Q (multiple  patches  allowed) . 

2.  The  exiting  panel  is  always  used  in  setting  up  the  partials 
associated  with  the  next  panel,  NQ. 

3.  A new  default  intersection  method  insures  an  appropriate 
intersection. 

4.  Once  an  intersection  point  is  found,  this  is  used  to  locate 
the  actual  next  panel,  NQ,  from  the  list  of  adjoining  panels. 

(a)  Points  in  any  interstices  require  that  each  possible 
next  peine  1,  NQ,  is  checked  for  streamline  success. 

(b)  The  new  NQ  must  have  peirticipated  in  the  calculation  in 
getting  across  peuiel  Q,  or  the  process  is  repeated. 

(c)  The  entering  point  is  required  to  be  on  the  edge  of  NQ. 

5.  A special  technique  is  used  for  elements  adjacent  to  a plane 
of  symmetry. 

6.  Streamlines  are  stopped  upon  approaching  a stagnation  region. 


Since  each  of  the  points  representing  a streamline  is  located  at  the 
boundary  between  two  quadrilaterals,  the  values  of  the  velocity  and  geo- 
desic curvatures  at  each  point  can  be  computed  separately  for  the  two 
quadrilaterals.  The  average  of  the  two  values  at  a point  is  taken  as 
the  value  at  that  point. 

The  streamline  is  first  traced  in  the  downstream  direction,  and  then  in 
the  upstream  direction.  After  the  entire  streamline  has  been  obtained, 
the  arc  length  along  the  streamline  is  computed  so  that  it  starts  at 
zero  at  the  upstream  end.  Points  which  are  very  close  together  are  com- 
bined. Such  points  occur  when  the  streamline  cuts  across  a corner  of 
a quadrilateral.  If  the  distance  between  two  points  is  less  than  one- 
eighth  the  distance  between  their  neighbor  in  front  and  their  neighbor 
in  back,  they  are  combined,  and  average  values  of  their  velocity,  posi- 
tion, and  curvatures  are  used  for  the  new  point. 
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SEARCH  ROUTINE 

Since  the  body  geometry  is  new  represented  by  a number  of  patches, 
each  containing  its  own  independent,  rectangular  network  of  panels,  it 
is  necessary  to  be  able  to  identify  panels  adjacent  to  each  other 
across  the  patch  edges;  the  user  describes  the  patches  and  the  particu- 
lar sides  of  the  patches  that  are  joined  together.  Further,  the  panel 
identifiers  along  any  side  of  a patch  are  calculated  internally.  In 
this  way,  searching  for  adjacent  panels  can  be  done  efficiently.  For 
example,  if  panels  adjacent  to  panel  N,  say,  on  side  2 of  patch  x are 
to  be  found,  then  only  panels  on  sides  of  the  patches  joined  to  side  2 
of  patch  x need  to  be  searched.  The  sketch  below  illustrates  the  me- 
thod used  to  select  adjacent  panels. 


To  test  whether  HQ  is  a neighbor  to  N,  the  following  quantities  sure  first 
formed  (see  sketch  for  vector  definitions) : 


(Mo  summation  is  implied  here  by  double  indices.)  For  NQ  to  be  a neigh- 
bor to  M, 


i * 1 or  2 


(38) 


(Da  < Ex  and  H.^  < E2) 

where  E^  and  E2  are  small  cutoff  numbers  which  are  discussed  later.  In 
the  sketch  above,  this  test  would  fail  for  i = 1 but  pass  for  i = 2,  and 
NQ  would  be  listed  as  a neighbor  to  N.  If  this  test  does  fail  for  i - 1 
and  i ■ 2,  then  a final  test  is  made  to  recover  NQ  panels  that  might 
completely  overlap  panel  N as  indicated  below: 


For  these  panels,  the  final  tests  that  must  be  satisfied  are 

(aL  * a2  < Ex  and  ^ < E2  and  H2  < E2>  (39) 


For  now,  and  E2  have  been  set  to  fairly  conservative  numbers  to  cover 
for  the  fact  that  neighboring  panel  edges  do  not  necessarily  lie  on  the 
same  line  because  panels  are  formed  in  the  mean  plane  of  the  four  corner 
points.  This  results  in  collecting  some  extra  panels  that  really  are 
not  adjacent  to  panel  N;  however,  the  subroutines  that  use  the  SEARCH 
information  have  automated  procedures  for  quickly  determining  the 
proper  adjacent  panel  needed  from  the  list  for  their  particular  problem 
(see  the  streamline  description,  for  example) . Actual  equations  for 
calculating  the  cutoffs  are 

Ei  - 1-1(PD»  * PV 


E 


2 


where  PD  is  the  projection  distance  used  in  generating  the  particular 
planar  panel. 
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HIND  TUNNEL 


The  capability  of  modeling  a body  with  winas  (at  zero  sideslip)  in  a 
wind  tunnel  has  been  added  to  the  program.  In  the  tunnel  option,  the 
tunnel  cross  section  is  defined  by  a set  of  comer  points.  In  between 
each  corner  point,  a set  of  lattice  points  will  be  generated  (see  sketch 
below) . 


3 Basic  Corner  Points 


Basic  Lattice  Points 


Thus,  the  user  inputs  each  basic  corner  point  and  the  number  of  lattice 
intervals,  NBT,  to  be  generated  by  linear  interpolation  between  this  cor- 
ner point  and  the  next.  The  default  value  is  1;  i.e.,  leaving  NBT  blank 
will  give  just  one  lattice  interval  between  this  comer  and  the  next  cor- 
ner input.  For  the  last  comer  point,  set  NB  to  a negative  integer;  i.e. 
-1. 

The  basic  lattice  points  (BLPY , BLPZ)  thus  generated  define  the  tunnel 
section  normal  to  the  free- stream  direction.  In  a set  of  numerical  cal- 
culations with  varying  a,  the  wind  tunnel  must  be  aligned  with  the  stream 
direction  each  time.  A pivot  point  (PX,  PZ)  must  be  clearly  defined  to 
keep  the  tunnel/model  relationship  correct.  The  PZ  value  must  have  the 
same  origin  as  the  defined  tunnel  section. 


Defined 

Section 
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In  rotating  tha  defined  section  normal  to  the  free  stream,  the  y-ordinate 
remains  unchanged.  A new  section  at  the  pivot  station,  located  at  the 
mid  station,  N,  is  then  generated  (so  as  not  to  lose  the  basic  one) . 
Sections  are  also  generated  at  the  leading,  L,  and  trailing,  T,  edges 
the  wind  tunnel  which  are  calculated  from  the  tunnel  half-length,  THL, 
input  by  the  user. 

The  streamwise  lattice  points  are  then  calculated  in  a cosine  spacing  to 
increase  lattice  density  at  the  pivot  station,  and  the  aerodynamic  in- 
fluence of  the  wind  tunnel  is  represented  by  a system  of  horseshoe  vor- 
tices as  indicated  in  the  following  sketch. 


With  wht,  defined  as  the  number  of  lattices  along  half  of  the  wind  tunnel 
then 


The  vortex  longitudinal  positions  are  at  equal  A0  intervals  with  the  con- 
trol points  displaced  A0/2.  Laterally,  the  control  points  are  placed 
midway  between  the  edges. 

The  tunnel  vortex  strengths  are  found  by  requiring  a sezo  normal  velocity 
at  each  tunnel  control  point.  Op  to  200  unknowns  can  be  used  to  represent 


the  wind  tunnel.  The  tunnel  influences,  along  with  all  cross-coupling 
terms,  have  been  assembled  into  the  basic  aerodynamic  matrix  as  indica- 
ted below. 
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ENGINE  INLET/EXHAUST  REPRESENTATION 


Because  the  engine  inlet  and  exhaust  of  a helicopter  configuration  are 
generally  an  integral  part  of  the  fuselage  flow  field,  it  is  important  to 
model  their  influence  in  the  analysis.  A model  that  has  been  used  suc- 
cessfully (Reference  6 ) is  shown  in  the  sketch  below. 


Slip  Surface 
(Vorticity  Sheet) 


Outflow  Plane 


Inflow  Plane 


6.  Geller,  E.W.,  Dvorak,  F.A.,  and  Risk,  M.,  "Calculation  of  the  Subsonic 
Flow  about  an  Advanced  Technology  Transport  and  a Modified  F-106B  Air- 
craft", Flow  Research  Note  No.  73,  Prepared  for  NASA-Lewis  Research 
Center,  Contract  NAS3-18341,  April  1975. 


The  inlet  and  exhaust  velocities  (mass  flow)  are  controlled  using  a pre- 
scribed normal  velocity  at  panel  control  points  in  the  inflow  and  outflow 
planes,  respectively.  The  exhaust  is  enclosed  within  a vortex  tube  that 
represents  the  slip  surface  between  the  exhaust  and  outer  flows.  A good 
approximation  to  the  vorticity  value  on  the  exhaust  tube  is  y = VQUt 

- V , which  assumes  that  the  exit  plane  pressure  is  that  of  the  free 
stream.  The  exhaust  model  is  therefore  basically  the  same  as  the  sepa- 
rated flow  model  in  this  program,  the  main  difference  being  the  sign  of 
the  vorticity,  since  the  former  encloses  a region  of  higher  energy  air 
and  the  latter  encloses  a region  of  lower  energy  air.  The  exhaust  case, 
however,  is  somewhat  easier  them  the  drag  case  because  the  vorticity  level 
is  known  approximately  at  the  outset,  being  related  to  the  engine  mass 
flow. 

A number  of  effluxes  are  allowed  in  the  program.  The  user  specifies  the 
points  around  the  edge  of  the  jet  outflow  plane,  the  length  of  the  efflux, 
and  the  vorticity  strengths  (i.e.,  VQOT  or  VJ£T,  approximately).  The 

program  generates  a set  of  streamwi.se  uniform  vorticity  panels  on  each 
efflux  tube  so  defined.  The  influence  of  each  efflux  tube  adds  directly 
to  the  right  side  of  the  boundary  condition  equations,  being  an  addi- 
tional onset  flow. 


APPLICATIONS 


A preliminary  version  of  the  program  DRAG  has  been  used  in  two  sepa- 
rate studies.  In  these  studies,  DRAG  has  been  exercised  in  a manner 
much  as  it  would  be  in  an  aircraft  development  program.  A preliminary 
version  of  the  code  was  used  which  had  all  of  the  features  of  the  final 
version  except  that  the  constant  source  body  singularities  were  retained. 
Execution  of  the  program  was  identical  to  the  version  with  the  higher- 
order  singularities,  since  these  will  be  completely  transparent  to  the 
user,  requiring  no  special  input.  All  the  features  of  the  program  were 
exercised,  cases  being  run  with  and  without  simulated  wind  tunnel  ef- 
fects, with  normal  flow  through  panels  to  simulate  inlet  and  exhaust 
flows,  with  simulated  efflux  plumes,  and  with  the  panel  addition  and 
subtraction  options  being  used.  In  addition,  although  not  required  for 
most  helicopter  configurations,  the  lifting  surface  portions  of  the  pro- 
gram have  been  thoroughly  evaluated.  All  of  this  work  involved  execu- 
tion of  the  geometry,  aerodynamic,  streamline,  and  boundary  layer  cal- 
culation elements  of  the  solution. 

The  program  application  was  carried  out  under  two  separate  studies. 

The  first,  funded  under  a subcontract  to  Hughes  Helicopters,  was  a study 
of  ways  to  reduce  the  drag  of  typical  utility  transport  and  armed  attack 
helicopters  through  a program  of  progressive  modification  of  the  hub  and 
pylon  region.  The  second  study  involved  the  NAVY/ARP A/Lockheed  X-Wing 
aircraft  and  was  funded  by  the  U.S.  Naval  Ship  Research  and  Development 
Center  (NSRDC) . In  both  studies,  the  analysis  was  being  used  to  analyze 
a baseline  configuration  and  then  to  explore  modifications  to  reduce 
the  drag. 

Panel  models  of  the  baseline  configurations  in  the  helicopter  hub/pylon 
drag  study  are  shown  in  Figures  3 and  4.  The  baseline  utility  trans- 
port helicopter  model  shown  in  Figure  3 was  paneled  by  Sikorsky  Air- 
craft under  a separate  Army-funded  study,  and  as  a consequence,  is  in 
the  Sikorsky  aerodynamic  program  format  (Ref.  7).  This  presents  no 
problem  since  the  latest  version  of  the  DRAG  program  has  the  flexibil- 
ity to  read  all  or  part  of  a model  in  the  Sikorsky  or  the  original 
WBAERO  input  format,  or  any  mix  of  the  two.  The  ability  to  read  in 
the  earlier  XYZ  program  format  has  been  retained  (Ref.  8).  The  attack 
helicopter,  shown  in  Figure  4,  was  set  up  from  the  beginning  in  the 
DRAG  program  input  format. 


7.  "Fuselage  Geometry  Definitions  Program",  Sikorsky  Aircraft  Report 
SER  50866,  September  1974. 

8.  Dawson,  C.W.,  and  Dean,  J.S.,  "The  XYZ  Potential  Flow  Program", 
Naval  Ship  Research  and  Development  Center  (NSRDC)  Report  3892, 
June  1972. 
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Using  the  flexibility  that  the  new  code  offers,  both  bodies  have  been 
set  up  using  multi-patch  paneling.  This  is  an  advance  over  the  earlier 
version  of  the  code  where,  in  order  to  track  streamlines,  the  paneling 
was  constrained  to  a rectangular  scheme  with  the  same  number  of  panels 
in  each  row  along  the  body.  Removal  of  this  constraint  means  that  com- 
plicated bodies  may  b=>.  paneled  efficiently.  An  example  of  a multi- 
patch body  with  the  patch  boundaries  outlined  is  given  in  Figure  1. 
Patching  bodies  gives  great  flexibility,  since  portions  such  as  the 
pylon  and  upper  body  may  be  changed  without  influencing  the  rest  of 
the  configuration.  Another  feature  of  the  study  was  that  inlet  and 
exhaust  mass  flows  were  to  be  simulated;  this  was  done  using  the  normal 
flow  option  in  the  porgram.  This  was  not  limited  to  engine  flows,  but 
was  also  used  to  model  the  flow  into  (and  out  of)  the  cavity  around  the 
rotor  shaft  and  to  simulate  escaping  cooling  air  from  the  attack  heli- 
copter avionics  compartment.  The  exhaust  plume  from  the  power  plants 
was  modeled  using  the  efflux  subroutine  now  available  in  the  program. 

An  important  part  of  the  helicopter  drag  study  was  the  testing  of  models 
at  each  stage  of  the  program  at  both  model  scale  (1/4)  and  almost  full 
(80%)  scale.  The  small  scale  testing  was  of  the  complete  aircraft, 
whereas  in  the  larger  scale  model  only  the  upper  fuselage,  pylon  and  hub 
were  used.  The  test  results  from  the  baseline  configurations  will  pro- 
vide data  for  correlation  of  the  analysis.  This  correlation  is  a 
necessary  prelude  to  the  application  of  the  analysis  in  the  later  stages 
of  the  contract  in  an  attempt  to  minimize  the  drag  of  both  configurations. 

Consequently,  the  analysis  was  used  not'only  to  model  the  complete  base- 
line aircraft  for  the  small  scale  tests,  but  also  to  model  the  partial 
model  on  the  ground  plane  for  the  large  scale  test.  Figure  5 shows  the 
partial  model  of  the  attack  helicopter  on  the  ground  plane.  The  wind 
tunnel  was  also  modeled  in  the  analysis,  but  is  not  shown  here  for 
clarity.  Also  evident  in  Figure  5 is  the  paneling  used  to  model  the 
rotor  hub  and  shaft.  While  a representative  model  of  the  hub  could  be 
constructed,  it  would  involve  detailed  paneling  and  separated  flow  model- 
ing beyond  the  scope  of  the  current  contract.  Since  its  presence  must 
be  modeled,  however,  it  was  done  in  the  following  manner.  From  earlier 
studies,  the  work  of  Sheehy  and  Clark  (Ref.  9)  is  typical,  it  has  been 
shown  that  a strong  correlation  exists  between  the  frontal  area  and  the 
drag  of  rotor  head  assemblies.  With  this  knowledge,  a crude  panel  model 
of  the  swept  volume  of  the  hub/shaft  may  be  constructed  and  normal  veloc- 
ities prescribed  which  produce  a momentum  deficit,  and  hence  a drag,  to 
match  the  experimental  data.  Used  in  conjunction  with  the  efflux  model 
attached  to  the  aft-facing  panels,  this  model  simulates  very  well  the 
presence  of  the  rotor  head  assembly  and  its  influence  on  the  surrounding 
surfaces. 


9.  Clark,  D.R.,  and  Sheehy,  T.W.,  "A  Method  for  Predicting  Helicopter  Hub 
Drag",  USAAMRDL-TR-75-48,  Eustis  Directorate,  U.S.  Army  A-ir  Mobility 
Research  and  Development  Laboratory,  Ft.  Eustis,  Va.,  January  1976. 
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Some  results  from  the  rotor  head  drag  study  are  given  in  Figures  6 
through  10.  In  Eigures  6 and  7,  the  centerline  pressure  distribution 
is  given  for  the  baseline  utility  and  attack  helicopters.  Very  clear 
in  both  plots  in  the  region  behind  the  rotor  head  and  down  the  aft 
pylon  is  a very  strong  unfavorable  pressure  gradient.  It  is  this  pres- 
sure gradient  which,  on  most  helicopter  configurations,  causes  the  mas- 
sive flow  separation  on  the  aft-facing  upper  surfaces.  The  unconnected 
points  in  both  plots  are  data  calculated  for  the  panels  on  the  hub  and 
shaft  model  and  on  the  inlets  around  the  shaft.  All  these  plots,  and 
the  geometry  plots  that  preceded  them,  are  produced  by  an  interactive 
data  review  package  that  was  developed  by  Analytical  Methods,  Inc., 
and  which  is  described  in  Reference  10.  Figure  8 shows  the  pressure 
distribution  around  a horizontal  cut  just  above  ground  plane  on  the 
partial  model  of  the  transport  helicopter  baseline.  The  distribution 
has  the  characteristic  shape  with  the,  again,  unfavorable  pressure 
gradient  as  the  flow  enters  the  transition  from  fuselage  to  tail  cone. 

Any  perturbation  in  this  region  of  unstable  flow  will  cause  separation 
and  drag.  Figure  9 is  a pressure  distribution  around  a body  station 
cut  close  to  the  crown  of  the  cockpit  canopy  on  the  utility  transport 
helicopter.  Of  particular  interest  is  the  plot  of  pressures  around  a 
buttline  cut  on  the  attack  helicopter  shown  in  Figure  10.  This  shows 
very  clearly  the  reduction  in  pressures  on  the  ground  plane  caused  by 
the  flow  accelerating  through  the  gap  between  ground  plane  and  engine 
nacelle.  The  acceleration  produces  effects  on  the  side  of  the  fuselage 
that  would  not  be  present  on  either  the  small,  full  model  or  on  the 
actual  aircraft.  Typical  of  what  will  happen  is  the  plot  of  the  stream- 
line trajectories  on  the  side  of  the  fuselage  shown  in  Figure  11.  This 
shows  how  the  streamlines,  numbered  4 and  5 in  the  illustration,  origi- 
nating on  the  front  fuselage  are  squeezed  down  by  the  nacelle  and  then, 
recovering,  shoot  up  onto  the  aft  pylon.  The  effect  of  this  flow,  in- 
teracting with  the  pylon  flow  in  a way  completely  unrepresentative  of 
the  real  world,  is  to  inhibit  the  aft  pylon  flow  separation.  Without 
the  prestudy  with  the  analysis  this  phenomenon  would  not  have  been 
discovered  without  a long  and  expensive  comparative  wind  tunnel  test. 

The  second  study  on  which  the  DRAG  program  has  been  exercised  is  the 
study  of  an  early  version  of  the  NAVY /ARP A/LOCKHEED  X-Wing  aircraft. 

The  aircraft  is  shown  without  wings  in  Figure  12,  and  in  the  full,  wings- 
fixed  configuration  in  Figure  13.  This  represents  perhaps  the  most 
challenging  case  that  has  been  executed  in  this  type  of  program  to  date. 
In  the  wings-off  case,  the  body  is  made  up  of  565  panels  in  a total  of 
26  different  patches.  This  large  number  of  patches  was  used  to  provide 
flexibility  in  changing  the  hub  fairing/pylon  configuration  to  reduce 


10.  Clark,  D.R. , "The  Use  of  Analytic  Tools  in  the  Design  and  Development 
of  Rotorcraft",  Presented  at  the  Fourth  European  Rotorcraft  and 
Powered  Lift  Aircraft  Forum,  Stresa,  Italy,  September  1978. 


Figure  9.  Pressure  Distribution  Around  a Body  Station  Cut  on  Transport 
Helicopter  Model. 


Figure  10.  Pressure  Distribution  Around  Buttline  Cut  on  Attack  Helicopter  Model 


STREAMLINE  DATA 


Streamline  Trajectories  on  Side  of  Attack  Helicopter  Fuselage 


the  drag.  A total  of  139  different  input  sections  were  used  to  build 
up  the  shape.  With  the  fully  lifting  wings  added,  the  number  of  panels 
is  increased  to  711. 
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The  object  of  the  study  is  to  analyze  the  flow  around  the  hub/pylon  for 
flight  Mach  numbers  up  to  0.7;  to  correlate  data  taken  in  an  earlier 
wind  tunnel  test;  and  to  suggest  detailed  changes  to  prevent  the  pos- 
sible occurrence  of  separated  flow,  and  hence,  drag.  Early  results  for 
the  wings-off  case  are  shown  in  Figure  14.  Here,  the  top  and  bottom 
centerline  pressures  are  plotted,  and  it  is  interesting  to  note  that 
even  with  a clean  rotor  head  fairing,  strong  unfavorable  pressure 
k gradients  exist  on  the  aft  pylon  which  will  likely  separate  the  local 

flow.  These  gradients  are  aggravated  in  the  lifting  case  where  the 
hub  fairing  suction  peak  coincides  with  the  aft  suction  peak  on  the 
trailing  wing.  A buttline  cut  through  the  wings  at  about  25%  span  is 
shown  in  the  next  figure.  The  strong  aft  loading  of  the  circulation- 
control  sections  used  in  this  aircraft  is  plainly  seen. 

The  exercise  of  the  program  in  the  two  applications  described  above  in 
parallel  with  its  final  development  has  provided  an  opportunity  to  ex- 
pose weaknesses  in  the  code  that  would  otherwise  not  have  emerged  until 
it  was  in  the  user's  hands.  While  the  basic  aerodynamics  code  has 
changed  little  from  the  earlier  version,  the  opportunity  to  model  com- 
plicated configurations  has  strengthened  the  modeling  capability  and 
has  led  to  wide  improvements  in  the  streamline  tracing  code. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  analysis  modified  to  permit  multi-patch  modeling  and  with  the  higher  order 
singularities  now  has  the  capability  to  predict  and  model  the  potential/viscous 
flow  around  generalized  configurations.  This  has  been  achieved  without  addi- 
tional complexity  in  either  input  or  execution.  The  program  has  been  exercised 
on  a number  of  configurations  representative  of  today's  technology,  and  no 
problems  have  emerged  from  running  these  test  cases. 

There  are,  however,  some  limitations  to  the  present  model,  removal  of  which 
would  further  enhance  the  program  capability.  The  first  limitation  is  that 
only  one  region  of  separated  flow  can  be  modeled.  This  may  be  acceptable  for 
some  shapes  with  a well  developed  base  flow  which  dominates  the  character- 
istics, but  for  detailed  modeling  on  some  types  of  helicopters  where  many 
regions  of  separated  flow  are  present,  further  program  development  is  required. 
The  second  limitation  lies  in  the  fact  that,  because  of  the  manner  in  which 
wing  paneling  is  generated,  information  required  to  calculate  the  streamline 
direction  is  not  available.  This  limits  calculation  of  viscous  effects  to  the 
body  panels.  For  helicopter  flow  field  modeling,  with  nonexistent  or  limited 
lifting  surfaces,  this  does  not  pose  a problem,  but  it  would  have  to  be  recti- 
fied if  calculation  of  the  separated  flow  on  the  wings  was  required. 

It  is  recommended  that  the  analysis  be  strengthened  by  inclusion  of  new 
elements  to  rectify  the  noted  restrictions. 
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APPENDIX  A 


PROGRAM  USER'S  GUIDE 

PROGRAM  INPUT  DESCRIPTION 

Card  Set  Input  Format 


1. 

Case  Title 

8A10 

2. 

RNB,  TRIPUP,  OPTION,  TRMAX,  XROUGH,  REFC, 

UIN,  XGEM 

8F10.0 

3. 

UKU,  DEN,  XO,  EPS,  XNROUGH 

7F10.0 

4. 

VRD(l)  ...  VRD (NROUGH) 

7F10.0 

5. 

XPRINT,  XSKIP , REFX,  REFZ,  CREF, 

PRINT,  CASE 

7F10.0 

6. 

NPANEL,  NBX,  MBX,  NWX,  MWX,  NVX, 
NS 

MVX,  I LX,  NPTS, 

1015 

7. 

TEXT  Card:-  GEOMETRY  INPUT 

8A10 

8. 

CASE,  PLOT,  SIM,  ISAVE,  PRINT 

5110 

9A. 

SINGPA,  NOPAN 

2110 

9B. 

X(I),  Y(I),  Z (I) 

7F10.0 

9C . 

NPAN  (1)  ...  NP AN (NOPAN) 

7110 

10A. 

NB 

110 

10S1. 

TYPE,  TITLE (I) 

18A4 

10S2. 

NOP,  NIBI 

2110 

10S3. 

XOF,  YOF,  ZOF,  XROT,  ZROT 

8F10.0 

10S4. 

AX (1)  , AY  (1)  , AZ (1)  ... 

...AX (4),  AY (4),  AZ(4) 

10X 

6F10.0 

10B. 

XBE,  YBE,  ZBE,  MB,  OPT,  FLAG 

3F10.0,  3110 

IOC. 

B (J)  , A(J)  , D(J) 

7F10.0 

10D. 

D(J),  B (J)  , A(J) 

3F12.0 

11A. 

NW,  KOORD 

2110 

54 


Cud  Set 

Input 

Format 

11B. 

XBE,  YBE,  ZBE,  CHRD,  ALF,  XAL,  MW,  OPT,  FLAG, 
DEL 

6F10, 

415 

11C. 

DELTA,  YO,  ZO 

7F10.0 

llO. 

B(J),  A(J) 

7F10.0 

11E. 

B(J),  D(J),  A (J) , C (J) 

5F12.0 

12. 

WAKE,  POINT 

7F10.0 

13A. 

XLP,  YLP,  ZLP 

7F10.0 

13B. 

XLP,  YLP,  ZLP 

7F10.0 

14A. 

NV 

110 

14B. 

XBE,  YBE,  ZBE,  MV,  OPT,  FLAG,  SIMPOT 

3F10.0,  4110 

14C. 

B (J)  , A(J)  , C(J),  D(J) 

7F10.0 

14D1. 

A,  B,  D 

7F10.0 

14D2. 

A,  B,  D 

7F10.0 

14D3. 

A,  B,  O 

7F10.0 

15. 

TEXT  Card:-  AERODYNAMIC  CALCULATIONS 

8A10 

16. 

NIT,  I EPS,  I TYPE 

3110 

17. 

COMPT,  SECT 

2110 

• 

00 

•H 

REFA,  REFL,  XOO,  X25 

7F10.0 

19. 

KUT,  NBV,  NV (1)  ...  NV (5) 

7110 

20. 

KOMPR,  POINTS 

2110 

2 1A. 

NHL,  PX,  PZ,  THL 

110,  3F10.0 

21B. 

YB,  ZB,  NBT 

2F10.0,  15 

22A. 

NEFF 

110 

22B. 

EFL(M) , EFGAM(M) 

7F10.5 

I 
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Card  Set 


Input 


Format 


I 


220. 

XE(I) , YE (I) , ZE(I) 

7F10.5 

23. 

NMA 

110 

24. 

MA 

7F10.0 

25. 

NAL 

110 

26. 

ALPHA,  BETA 

7F10.0 

27. 

IS 

110 

28. 

IA,  IE 

2110 

29. 

IREI 

110 

30. 

11(1),  11(2),  ...  ETC. 

1415 

31. 

DELY,  REFL,  XLE 

7F10.0 

32. 

TEXT  CARD:-  STREAMLINE  CALCULATIONS 

8A10 

33. 

ENDFILE 

110 

34. 

NLIN,  NSP(l),  ...  NSP(NLIN) 

1315 

35A. 

NJT 

1315 

35B. 

JOINT (I, J) 

1315 

36. 

6/7/8/9. 

56 
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Description  of  Input  Variables 


Card  1:-  General  Identification, 
information  in  Columns  1-80. 

Card  2:-  General  Parameters. 

Column  Variable  Value 

1-10  RNB  arbitrary 


11-20 


21-30 


TRIPUP 


OPTION 


31-40 


41-50 


TRMAX  arbitrary 


XROUGH 


51-60 


REFC  arbitrary 


Card  1 contains  any  desired  identifying 


Description 

Reynolds  number  based  on  reference 
chord  and  free-stream  velocity 
U^C/V  x (10'6) 

N.B.  If  the  input  geometry  is  in  di- 
mensional units,  the  Reynolds  number 
will  be  input  per  unit;  i.e.,  Re/inch 

Trip  location  (x/c) 

NO  tripping  desired 

Boundary  layer  will  trip  where  speci- 
fied by  TRIPUP 

Program  will  test  on  boundary  layer 
momentum  thickness  at  trip  location. 

If  Rq  < 200,  trip  location  will  be  re- 
positioned to  point  where  Rq  >_  200. 

This  deters  the  user  from  specifying 
an  unrealistically  early  trip  location. 

Maximum  number  of  iterations  between  in- 
viscid  and  boundary  layer  modules 

Standard  boundary  layer  calculation  for 
smooth  surfaces 

Boundary  layer  calculation  with  area 
suction  (needs  Cards  3 and  4) 

Boundary  layer  calculation  with  surface 
roughness  (needs  Cards  3 and  4) 

Reference  chord  in  inches  for  determina- 
tion of  surface  distance  in  INSPAN 


Variable 


Value 


Column 


Description 


61-70 


(JIN 


arbitrary 


Free-stream  velocity  in  feet  per  second 
for  use  in  INSPAN 


71-80 


XGEM 


0. 


Standard  case;  program  will  execute 
for  complete  case 


1.  Program  will  assemble  and  print  out 

geometry  variables  only  and  stop 

N.B.  REFC,  UIN  and  XGEM  currently  not  available. 


Card  3:-  Roughness  Parameters. 


Column 

Variable 

Value 

Description 

1-10 

UKU 

30.-70. 

Value  of  roughness  Reynolds  number 
where  effect  of  roughness  is  indepen- 
dent of  viscosity  (see  Ref.  1) 

11-20 

DEN 

2.-100. 

Density  of  roughness  element  spacing 
(For  sand  grain  roughness,  use  3.175) 

21-30 

XO 

.005-. 02 

Initial  value  of  skin  friction  coeffi- 
cient (.007  is  usual  starting  value) 

31-40 

EPS 

.0001 

Error  bound  in  skin  friction  calcula- 
tion 

41-50 

XNRQUGH 

arbitrary 

Number  of  input  values  of  suction  velo- 
city, v ^/U,  or  roughness  heights,  k/c 

Card  4s- 

Roughness  or 

Suction  Distribution. 

Column 

Variable 

Value 

Description 

1-10 

• 

VRD(l) 

• 

arbitrary 

(floating 

point) 

Roughness  height  or  area  suction  velo- 
city (k/c  or  V/Vm) 

61-70 


VRD (NROUGH) 


N 


•a 


Card  5:—  General  Parameters 


Column 

Variable 

Value 

Description 

1-10 

XPRINT 

. 0. 

Suppresses  printing  of  cross-flow  inte- 
gral thickness  from  IBL 

1. 

Extra  printing 

11-20 

XSKIP 

1. 

Print  option  in  INTEGRAL;  every  integra- 
tion step  is  printed 

10. 

Every  tenth  step  is  printed 

21-30 

REFX 

arbitrary 

Reference  (x/c)  location  for  calculation 
of  moment  coefficient 

31-40 

REPZ 

n 

Reference  (z/c)  location  for  calculation 
of  moment  coefficient 

41-50 

CREF 

l. 

Redundant 

51-60 

PRINT 

0. 

II 

61-70 

CASE 

l. 

Number  of  angle-of-attack  or  Mach  number 
variations  for  a given  geometry.  Cur- 
rently limited  to  1 

Card  6s- 

Variable  Dimension  Options 

Column 

Variable 

Value 

Description 

1-5 

NPANEL 

.LE.1500 

Number  of  source  panels  and  vortex  lat- 
tices on  one  side  of  plane  of  symmetry 

6-10 

NBX 

2-70 

Number  of  body  sections 

11-15 

MBX 

3-60 

Maximum  number  of  input  points  for  any 
body  section 

16-20 

NWX 

2-40 

Number  of  wing  sections 

21-25 

MWX 

3-59 

Maximum  number  of  ordinates  at  any  wing 
section 

26-30 

NVX 

0-40 

Number  of  streamwise  vortices  in  body 
lattice 

31-35 

MVX 

2-60 

Number  of  bound  vortices  in  body  lattice 

59 


Column 


Variable 


Value 


] 


Description 


36-40 

ILX 

2-35 

Sum  of  wing  and  body  lattices 

41-45 

NPTS 

.LE.1500 

Nusber  of  off-body  points 

46-50 

NS 

.LE.8 

Number  of  streamlines 

N.B. 

Great  care  should  be  taken  in  preparing  this  card  since  experience 
has  shown  that  errors  here  are  the  most  common  cause  of  aborted 
runs.  Too  large  a number  selected  causes  machine  capacity  problems; 

too  small  gives  inadequate  array  space  and  overwrites. 
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Geometry  Input  Cards If  the  configuration  is  symmetrical  about  the  x,z 
plane,  geometrical  input  is  required  for  only  one  side  of  the  configura- 
tion. The  convention  used  herein  is  to  present  that  half  of  the  configu- 
ration lying  on  the  positive  y side  of  the  x,z  plane.  If  the  configura- 
tion is  not  symmetric,  complete  geometrical  input  is  required. 


Card  7i-  General  Identification.  Card  7 contains  any  desired  identifying 
information  in  Columns  1-80. 


Card  61-  Options  Card. 

Coliann  Variable  Value 

10  CASE  1 

2 


Description 
Isolated  body  only 
Isolated  wing  only 


3 

20  PLOT  0 

1 

30  SIM  0 

1 


Wing-body  combinations 

No  plot  output  (currently  this  is  the 
only  option  available) 

Plot  output  requested 

Configuration  symmetric  about  x,z  plane; 
panel  geometry  required  on  one  side  only 
(normal  case) 

Configuration  synaetric  about  x,z  plane. 
Panel  geometry  input  required  on  one 
side  only;  panel  geometry  output  calcu- 
lated for  both  sides.  (Used  when  analy- 
zing synaetric  configuration  in  yaw.) 

Not  currently  available 


J 

1 

•i 


Column 

Variable 

Value 

Description 

30 

SIM 

-1 

Asymmetric  configuration.  Panel  geo- 
metry input  required  for  both  sides 

40 

I SAVE 

0 

Geometry  and  influence  coefficient 
matrices  not  saved 

1 

Geometry  and  influence  coefficient 
matrices  saved  in  previous  run  to  be 
used  (TAPES  11  and  18  must  be  requested) 

-1 

Geometry  and  influence  coefficient  ma- 
trices to  be  saved  on  TAPES  11  and  18 

50 

PRINT 

0 

Normal  output 

1 Optional  output  1.  Includes  panel  geo- 
metry, coordinate  transformation  matri- 
ces, and  panel  forces  and  moments 

2 Optional  output  2.  Panel  velocity  com- 
ponents and  influence  coefficients. 
Requires  large  line  count  limit 

3 Optional  output  3.  Hie  aerodynamic 
influence  coefficient  matrix,  the  right 
side  of  the  matrix  equation,  and  all 
solution  iterations 

4 Optional  output  4.  This  option  prints 
out  the  successive  solution  iterations 
only 

N.B.  The  normal  output  is  always  printed  in  addition  to  any  optional  out- 
put selected. 


Card  9A:-  Single  Panel  Control  Card. 
Column  Variable  Value 

1-10  SINGPA  0 N 


Description 


No  single  panel  input;  omit  Card  Set  9B, 
continue  reading  input  cards.  Corner 
point  coordinates  of  this  panel  follow 
on  Card  Set  9B 


11-20  NOPAN  arbitrary  Number  of  panels  to  be  deleted.  If  non- 

integer zero,  panel  indices  follow  on  Card  Set  9C 


0 


Card  9B»-  Panel  Corner  Point  Input. 


Column 

Variable 

Value 

Description 

1-10 

X(X) 

arbitrary 

(floating 

x-coordinate  of  corner  I 

point) 

11-20 

YU) 

N 

y-coordinate  of  corner  I 

21-30 

XU) 

M 

z-coordinate  of  corner  I 

Repeat  Card  9B  four 

times,  once 

for  each  corner  of  the  panel 

Card  Set  9C:-  Indices  of  Deleted  Panels.  NOPAN  indices  of  deleted  panels 
are  read  (7110  format)  if  NOPAN  > 0 on  Card  9A.  A maximum  of  100  panels 
may  be  deleted.  Ming  and  body  vortex-lattice  control  panels  may  not  be 

deleted. 

HARMING:  Great  care  should  be  taken  when  exercising  this  option  when 

going  to  a streamline  trace  since  dislocation  of  the  row,  column 
and  panel  counting  can  occur.  Its  use  is  not  recommended  until 
the  user  is  completely  familiar  with  program  operation. 


Card  Set  10:-  Body  Panel  Input.  This  card  set  allows  the  body  panels  to 
be  calculated  automatically  from  the  section  geometry  data.  Five  options 
are  available  for  inputting  the  section  geometry.  The  XYZ  program  input 
referred  to  below  conforms  with  the  format  of  Reference  2.  The  Sikorsky 
input  conforms  to  that  outlined  in  Reference  9 for  input  to  program  Y 179C 
with  minor  additions  outlined  below  (omit  this  Card  Set  if  CASE  = 2 on 
Card  8) . 

Card  lOAt-  Number  of  Body  Sections. 

Column  Variable  Value  Description 

1-10  NB  arbitrary  Number  of  body  sections  (2  <.  NB  140) 

integer  A negative  number  indicates  Sikorsky 
input  format 


Card  Set  lOSi-  Alternate  Sikorsky  input.  This  card  set  is  entered  by 
loading  NBflQA)  with  a negative  number.  The  value  is  not  significant. 
Standard  WBAERD/DRAG  input  may  be  resumed  after  reading  Sikorksy  input 
by  reloading  NB(IOA)  with  the  appropriate  positive  value.  Panels  should 
be  grouped  for  input  in  regular  arrays,  as  in  a DRAG  program  patch,  each 
patch  of  panels  being  separated  by  a Type/Title  card  (10S1),  a panel  count 
card  410S2)  and  a local  origin  and  patch  rotation  card  (10S3) . To  ter- 
minate a load  in  the  Sikorsky  input  format,  simply  repeat  the  10S  Type 
Card  Set  appropriately. 


Card  10S1 : - Type/Title. 

Column  Variable  Value 

1-4  TYPE  NL,S 


BLANK 

END 

TITLE  alphabetic 


Description 

Program  expects  nonlifting  (body) 
panels,  symmetric  about  centerline. 
(This  is  currently  the  only  option 
available) . Use  on  first  input  set 
only 

Normal  input  sets 
Terminates  Sikorsky  input 
Any  title 


Card  Set  10S2:-  Panel/ 


Countir 


Column 

Variable 

Value 

Description 

1-10 

NOP 

arbitrary 

integer 

Number  of  panels  in  patch 

11-20 

MBI 

arbitrary 

integer 

Number  of  panels  in  each  column  of 
patch 

Card  10S3: 

- Offsets 

and  Rotation. 

Column 

Variable 

Value 

Description 

1-10 

XOF 

arbitrary 

(floating 

point) 

x-coordinate  of  patch  local  origin 

11-20 

YOF 

•• 

y-coordinate  of  patch  local  origin 

21-30 

ZOF 

m 

z-coordinate  of  patch  local  origin 

I 


Column 


Variable 


Value 


r 


r 


31-40 


41-50 


Description 


XROT 


ZROT 


arbitrary  Angle  in  degrees  which  input  points 
(floating  should  be  rotated  about  x-axis 
point) 


Angle  in  degrees  which  input  points 
should  be  rotated  about  z-axis 


Card  Set  10S4:-  Panel  Corner  Points. 

Card  10S4.1. 


Column 

Variable 

Value 

Description 

1-10 

BLANK 

May  be  used  for  panel  identification 

11-20 

AX(1) 

arbitrary 

(floating 

point) 

Ordinates  of  first  corner  point  on  first 
panel 

21-30 

AY  (1) 

«• 

M 

31-40 

AZ(1) 

N 

M 

41-50 

AX  (2) 

n 

Ordinates  of  second  corner  point  on 
first  panel 

51-60 

AY  (2) 

m 

It 

61-70 

AZ(2) 

M 

N 

I 

i 

i 

t 

j 


Card  10S4 
Column 
1-10 
11-20 

21-30 

31-40 


Variable 
BLANK 
AX  (3) 

AY  (3) 
AZ(3) 


Value 


arbitrary 

(floating 

point) 

N 


Description 


Ordinates  of  third  corner  point  on  first 
panel 

n 


Ordinates  of  fourth  corner  point  on 
first  panel 


II 

I 
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41-50 


AX  (4) 


r 


Column 


Variable 


Value 


Description 


51-60 


61-70 


AY  (4) 


AZ(4) 


arbitrary 

(floating 

point) 


Ordinates  of  fourth  corner  point  on 
first  panel 


Repeat  Card  Set  10S4  NOP  times  to  cctnplete  the  input  for  each  patch. 
Repeat  Card  Set  10S  for  each  patch  to  be  loaded. 

To  terminate  Sikorsky  input,  load  Card  10S1  with  TYPE  = END. 


Card  10B:- 


Section  Geomet 


Column 


Variable  Value 


Description 


11-20 


21-30 


31-40 


arbitrary  x-coordinate  of  origin  of  body  section 
(floating  coordinate  system,  except  blank  when 
point)  XYZ  program  input  format  is  used 

" Similarly  the  y-coordinate 

" Similarly  the  z-coordlnate 

arbitrary  Number  of  input  points  on  section 

integer  (3  £ MB  £ 60) 

If  MB  < 0,  XYZ  program  input  format  re- 
quested 

0 Body  section  geometry  input  by  y-z  co- 
ordinates on  Card  Set  10C  and  10D 

1 Body  section  geometry  same  as  preceding 
body  section  - Card  Set  10C  or  10D 
emitted.  Note:  YBE  and  ZBE  are  addi- 
tive to  preceding  values 

2 Body  section  geometry  input  in  polar 
coordinates,  r,  6,  on  Card  Set  10C 

3 Body  of  revolution,  section  geometry 
input  as  section  radius  and  theta  incre- 
ment on  Card  Set  10C 


* 


Column 

Variable 

Value 

Description 

60 

FLAG 

0 

Normal  body  section 

1 

Terminal  body  section  (end  of  current 
body  panel  network) 

2 

End  of  rectangular  grid 

Card  Set 

IOC:-  Body  Section  Coordinates  (Normal  Input) . 

Column 

Variable 

Value 

Description 

1-10 

B(J) 

arbitrary 

(floating 

point) 

y-coordinate  of  point  J if  OPT  = Of  or 
angular  coordinate  (in  degrees)  of  J 
if  OPT  = 2;  or  increment  angle,  A 8,  in 
degrees  if  OPT  = 3 

♦ 

11-20 

A(J) 

N 

z- coordinate  of  point  J if  OPT  = 0;  or 
r-coordinate  of  point  J if  OPT  = 2;  or 
body  section  radius  if  OPT  = 3 

21-30 

D(J) 

n 

A x shift  of  point  J if  OPT  = 0 or  2 

Card  Set 
OPT  =*  3; 

IOC  contains  MB  cards  if  OPT  = 0 or  2;  contains  only  1 card  if 
and  is  omitted  if  OPT  * 1 or  MB  < 0 on  Card  10B. 

Card  Set 

100:-  Alternate  ZYZ  Input 

Column 

Variable 

Value 

Description 

1-12 

D(J) 

arbitrary 

(floating 

point) 

x-coordinate  of  point  J 

13-24 

B(J) 

If 

y-coordinate  of  Point  J 

25-36 

A(J) 

If 

z-coordinate  of  point  J 

This  card  set  is  omitted  unless  MB  < 0 in  card  10B. 

Note:  Repeat  Card  10B  and  Card  Sets  IOC  or  10D  NB  times  to  complete  Card 
Set  10. 
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Card  11A;-  Number  of  Wing  Sections. 


Column 

Variable 

Value 

Description 

1-10 

NW 

arbitrary 

integer 

Number  of  wing  sections  (2  <_  NW  <_  40) 

20 

KOORD 

1 

Wing  section  ordinates  input  in  percent 
of  local  chord 

2 

Wing  section  ordinates  input  are  not 
normalized 

Card  11B:- 

Wing  Section  Geometry. 

Column 

Variable 

Value 

Description 

1-10 

XBE 

arbitrary 

(floating 

point) 

x-coordinate  of  origin  of  wing  section 
coordinate  system 

11-20 

YBE 

IV 

Similarly  the  y-coordinate 

21-30 

ZBE 

If 

Similarly  the  z-coordinate 

31-40 

CHRD 

II 

Chord  length  section 

41-50 

ALF 

arbitrary 

(floating 

point) 

Section  twist  angle  (degrees)  or  flap 
rotation  angle.  (Degrees  positive  for 
positive  flap  deflection) 

51-60 

XAL 

•I 

Center  of  twist  or  rotation  in  percent 
chord 

61-65 

MW 

arbitrary 

integer 

Number  of  coordinates  in  section 
(5  < MW  < 59) 

70 

OPT 

0 

Wing  section  ordinates  to  be  used  from 
card  set  11D 

1 

Wing  section  ordinates  same  as  preced- 
ing section  - card  set  11D  omitted 

75 

FLAG 

0 

Normal  case  - surface  vorticity  calcula- 
ted automatically 

1 

Terminal  wing  section  (end  of  current 
wing  panel  network) 

Column 


Variable 


Value 


75  FLAG 


80  DEL 


Description 

2 No  vortex-lattice  panels  calculated 

for  this  section 

3 The  coordinates  of  the  last  bound  vor- 
tex in  the  vortex  lattice  are  read  in 
on  Card  13  for  this  section 

0 No  wing  dihedral 

1 Dihedral  input  on  Card  Set  11C 


Card  Set  11C;-  Wing  Dihedral  Input 
Column  Variable  Value 


1-10  DELTA  arbitrary 

(floating 

point) 


Description 


Dihedral  angle  (degrees) 


11-20  VO 


y-  and  z-ordinates  of  axis  of  rotation 
of  wing  panel 


21-30  ZO 

N.B.  Omit  Card  Set  11C  if  DEL  = 0 in  Card  11B. 


Card  Set  IIP:-  Wing  Section  Coordinates. 

Column  Variable  Value  Description 

1-10  B (J)  arbitrary  x-coordinates  of  point  J 

(floating 

point) 

11-20  A(J)  " z-coordinate  of  point  J 

N.B.  Card  Set  11D  contains  MW  cards  of  OPT  = 0,  and  is  omitted  if  OPT  = 
Repeat  Card  11B  and  Card  Sets  11C  and  11D  NW  times  to  complete  Card 
Set  11. 
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Card  Set  12;-  Vortex-Lattice  Control  Point  Location. 


Column 


Variable 


Value 


Description 


1-10 


WAKE 


arbitrary 

(floating 

point) 


Extension  of  vortex  lattice  into  wake 
in  percent  chord  (usually  100.) 


11-20 


POINT 


Location  of  vortex-lattice  control  point 
in  percent  chord  behind  trailing  edge 
(usually  .1) 


N.B.  These  values  are  not  used  if  FLAG  - 3 on  Card  11B. 


Card  Set  13:-  Relocation  of  Vortex-Lattice  Terminal  Points.  This  Card  Set 


is  omitted  unless  FLAG  - 3 on  Card  11B.  For  each  wing  section  having  FLAG 
- 3,  two  additional  cards  are  required  to  specify  the  terminal  points  of 
the  streamwise  vortices. 


Card  Set  13A:-  Inboard  Terminal  Points. 


Column 

Variable 

Value 

Description 

1-10 

XLP 

arbitrary 

(floating 

point) 

x-coordinate  of  inboard  edge  of  lattice 
terminal  point 

11-20 

YLP 

N 

y-coordinale  of  inboard  edge  of  lattice 
terminal  point 

21-30 

ZLP 

M 

z-coordinate  of  inboard  edge  of  lattice 
terminal  point 

Card  Set  13B:-  Outboard  Terminal 
edge  of  lattice  terminal  point. 

Points.  Same  as  Card  13A  for  outboard 

j t 

Card  Set 

14:-  Body  Vortex-Lattice  Input.  This  Card  Set  allows  additional 

vortex  lattices  to  be  located  inside  the  body  of  wing-body  combinations, 
and  is  omitted  if  CASE  < 3 on  Card  8. 
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Network 


Column  Variable  Value  Description 

1-10  NV  arbitrary  Number  of  streamwise  vortices  in  body 

integer  vortex- lattice  network  (NV  < 40) 

• 1 

Note:  The  sum  of  all  wing  and  body  vortex  lattices  may  not  exceed  35. 


I 


Card  Set  14B:-  Vortex-Lattice  Geometry. 


Column 

Variable 

Value 

Description 

- 

1-10 

XBE 

arbitrary 

(floating 

point) 

x-coordinate  of  origin  of  streamwise 
vortex 

11-20 

YBE 

n 

y- coordinate  of  origin  of  streamwise 
vortex 

21-30 

ZBE 

n 

z-coordinate  of  origin  of  streamwise 
vortex 

31-40 

MV 

arbitrary 

integer 

Number  of  bound  vortices  in  lattice 

2 < MV  < 60 

41-50 

OPT 

0 

Vortex-lattice  points  to  be  read  from 

Card  Set  14C 

1 

Vortex- lattice  points  same  as  preceding. 
Omit  Card  Set  14C 

2 

Optional  vortex-lattice  control  panel 
coordinates  read  on  Card  14D3 

51-60 

FLAG 

0 

Normal  case  - vortex-lattice  panels 
calculated 

1 

Terminal  vortex  of  current  body  vortex- 
lattice  network 

» 

2 

Corner  points  of  control  point  panel  to 
be  read  on  Cards  14D2  and  14D3  (used 
when  arbitrary  control  point  is  desired) 

* 

Value 


Synnetry  option  specified  on  Card  8 
enforced  for  this  vortex 


Synnetry  option  ignored  for  this  vor- 
tex lattice  (used  for  inserting  vortex- 
lattice  networks  in  vertical  tails  lo- 
cated in  x,z  plane) 


ordinates . Two  or 
minal  point  of  the 
control  point  panel 

three  additional  cards  jure  required  to  specify  the  ter- 
streamwise  vortex,  and  the  comer  points  of  the  lattice 

Card  14D1. 

Column  Variable 

Value 

Description 

1-10  B 

arbitrary 

(floating 

point) 

x-coordinate  of  terminal  point  of  stream- 
wise  vortex 

11-20  A 

M 

z-coordinate  of  terminal  point  of  stream- 
wise  vortex 

21-30  D 

n 

Ay  shift  of  terminal  point  of  streamwise 
vortex 

1-10 

B(J) 

arbitrary 

(floating 

x-coordinate  of  point  J 

point) 

11-20 

A(J) 

N 

z-coordinate  of  point  J 

21-30 

C(J) 

It 

Vortex- lattice  strength  at  point  J 

31-40 

D(J) 

n 

Ay  shift  of  point  J 

Card  Set  14C 

contains  MV  cards  if 

OPT  *0,  and  is  omitted  if  OPT  = 1 

on 

Card  14B. 

Control  Set 

14D:- 

Vortex-Lattice 

Terminal  Point  and  Control  Point  Co- 

Note:  This  point  also  defines  the  upstream  corner  of  the  control  point 
panel  if  FLAG  ft  2 on  Card  14B. 


Card  14d2 : - Same  as  Card  14D1,  containing  the  coordinates  of  the  down- 
stream corner  of  the  control  point  panel  if  FLAG  ft  2 on  Card  14B.  If 
FLAG  * 2 on  Card  14B,  this  card  contains  the  coordinates  of  the  upstream 
corner  of  the  control  point  panel. 


Card  14P3:~  If  FLAG  = 2 on  Card  14B,  this  card  contains  the  coordinates 
of  the  downstream  corner  of  the  control  point  panel  in  the  same  format  as 
Card  14D1.  Omit  this  card  if  FLAG  ft  2 in  Card  14B1. 

Note:  Repeat  Card  14B  and  Card  Sets  14c  and  14D  NV  times  to  complete  Card 
Set  14. 

Aerodynamic  Input  Cards:-  The  configuration  panel  geometry  is  transferred 
to  the  aerodynamic  section  of  the  program  by  TAPE  11.  Additional  aerody- 
namic input  cards  required  are  described  below. 


Card  15:-  Case  Identification  Card.  Card  15  contains  any  der  x J ise 
identification  in  Columns  1-80. 


Card  16:-  Iteration  Option  Card. 
Column  Variable  Value 


Description 


NIT 

arbitrary 

integer 

Maximum  number  of  iterations  (15  - 25) 

IEPS 

N 

Exponent  of  10  setting  limiting  value 
for  residue  of  iterative  solution  (-3 
or  -4  recommended) 

ITTPE 

1 

Blocked  Jacobi  iteration  procedure 

2 

Blocked  Gauss-Seidel  iteration  procedure 

3 

Blocked  Gauss-Seidel  with  controlled 
successive  over-relaxation 

4 

Blocked  Gauss-Seidel  with  successive 
over-relaxation 

Option  4 has  not  been  completely  checked  out,  but  will  be  made  avail- 
able at  a later  date. 
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Card  17  >-  Configuration  Options. 
Column  Variable  Value 

0 


Description 


1-10 


COMPT 


Forces  and  moments  calculated  for  com- 
plete configuration 


arbitrary 

integer 


Forces  and  moments  calculated  on  com- 
ponents. Panel  indices  of  each  com- 
ponent follow  on  Card  26 


11-20 


SECT 


0 

1 


No  wing  section  forces  and  moments 


Wing  section  forces  and  moments  calcu- 
lated. Wing  section  indices  follow  on 
Card  25,  panel  indices  in  each  section 
on  Card  26,  and  section  reference 
lengths  on  Card  29 


Forces  and  moments  calculated  on  sub- 
sections. The  number  of  subsections 
follow  on  Card  25,  the  number  of  panel 
groi$>s  on  Card  27,  the  panel  indices 
in  each  group  on  Card  28,  and  subsection 
reference  lengths  on  Card  29 


Card  18:- 

Reference 

Parameters . 

Column 

Variable 

Value 

Description 

• 

1-10 

REFA 

arbitrary 

(fixed 

point) 

Reference  data 

1 

11-20 

REFL 

n 

Reference  chord  (MAC) 

1 

21-30 

XOO 

it 

Axial  distance  of  leading  edge  of  MAC 
from  origin 

j 

l 

31-40 

X25 

it 

Axial  distance  of  quarter  chord  of  MAC 

from  origin 

Default  option — all  variables  set  to  1.  internally  (leave  card  blank). 


i 
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Card  19 t- 

Configuration  Lift  Option » 

a 

Column 

Variable 

Value 

Description 

• 

1-10 

KUT 

0 

Nonlifting  configuration,  no  vortex-lat- 
tice Kutta  condition  imposed 

1 

Lifting  configuration,  vortex- lattice 

Kutta  condition  imposed 

-1 

Wing  vortex  lattice  extends  through  body 
having  same  strength  as  adjacent  wing 
vortex  lattice 

• 

j 

11-20 

NBV 

arbitrary 

integer 

Number  of  body  vortices  (NBV  £ 5) 

21-30 

NV(1) 

n 

Number  of  wing  vortices  associated  with 
body  vortex  1 

31-40 

• 

NV(2) 

• 

N 

Number  of  wing  vortices  associated  with 
body  vortex  2 

• 

61-70 

• 

NV(5) 

n 

Number  of  wing  vortices  associated  with 
body  vortex  5 

Card  20:- 

Compressibility  Rule  Option. 

Column 

Variable 

Value 

Description 

1-10 

KOKPR 

1 

Gothert  Rule  1 selected 

2 

Gothert  Rule  2 selected 

l 

11-20 

POINTS 

0 

One  body  component  case,  no  off-body 

points 

1 Flap  case,  off-body  points  will  be 

calculated 
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Card  Set  21:-  Wind  Tunnel  Input.  The  manner  in  which  the  wind  tunnel  is 
modelled  is  outlined  in  the  main  body  of  the  report. 


Card  Set  21B:-  Tunnel  Cross-Section. 

Column  Variable  Value  Description 

1-10  YB  arbitrary  y coordinate  of  point  defining  an  element 

(floating  of  tunnel  cross-section 

point) 

11-20  ZB  " 2 coordinate  of  point  defining  an  element 

of  tunnel  cross-section 

H 

21-25  NBT  arbitrary  Number  of  panels  into  which  line  between 

integer  this  point  and  next  point  is  cut 

-1  Last  point  defining  section 


Repeat  Card  Set  21B  as  many  times  as  there  are  straight  line  elements  in 
tunnel  cross-section. 


Card  22A:-  Number  of  Effluxes 


Column 

Variable 

Value 

Description 

1-10 

NEFF 

arbitrary 

integer 

(<5) 

Number  of  separate  effluxes 

Card  22B: 

- Efflux  Length  and  Strength. 

Column 

Variable 

Value 

Description 

1-10 

EFL(  ) 

arbitrary 

(floating 

point) 

Efflux  length 

11-20 

EFGAMt  ) 

M 

Efflux  strength 

Card  22C: 

- Efflux  Cross-Section. 

Column 

Variable 

Value 

Description 

1-5 

NPAN 

arbitrary 

integer 

Number  of  body  panel  edges  defining 
efflux 

(NPAN*NEFF<60) 

Card  Set 

220:-  Efflux  Cross-Section  Details. 

Column 

Variable 

Value 

Description 

1-10 

XE(1) 

arbitrary 

(floating 

point) 

ordinates  of  points  defining  efflux 
opening 

11-20 

YE  (I) 

n 

If 

21-30 

ZE(1) 

n 

M 

Repeat  Card  Set  22D  NPAN+1  times  Cor  each  efflux.  Repeat  Cards  22B  and 

22C  and  Card  Set  220  for  each  efflux. 

Card  23:-  Number  of  Mach  Nuabers. 

Column 

Variable 

Value 

Description 

1-10 

NMA 

arbitrary 

integer 

Number  of  Mach  numbers  following  on  Card 
Set  24 
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Card  29:-  Number  of  Panels  in  Subsections 


Column 

Variable 

Value 

Description 

1-10 

IREI 

arbitrary 

integer 

Number  of  panels  in  subsections.  Omit 
if  SECT  < 2 on  Card  17 

Card  30:- 

Subsection 

Panel  Indices. 

Column 

Variable 

Value 

Description 

1-5 

II  (1) 

arbitrary 

integer 

Panel  indices  of  all  panels  in  subsec- 
tion; omit  if  SECT  < 2 on  Card  17 

6-10 

11(2) 

•« 

11-15 

11(3)  ...i 

etc. 

Card  31:- 

Reference  Lengths. 

Column 

Variable 

Value 

Description 

1-10 

DELY 

arbitrary 

integer 

Width  of  section 

11-20 

REFL 

II 

Reference  length  of  section 

21-30 

XLE 

•9 

Moment  reference  point  of  section 

N.B.  Cards  27  - 31  must  be  repeated  for  each  angle-of-attack  or  yaw,  if 
section  data  requested. 


Streamline  Input  Cards:-  Each  streamline  selected  is  constrained  to  pass 
through  a chosen  panel.  This  avoids  the  possibility  of  the  selected 
streamlines  bunching  together  in  the  separated  flow  region. 


Card  32:-  Case  Identification  Card.  Card  32  contains  identification  for 
the  streamline  calculations. 


Card  33:-  Streamline  Option  Card. 
Column  Variable  Value 

1-10  ENDFILE  -1 


Description 

Program  switches  out  of  WBAERO  onto  the 
streamline  overlay 
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Card  34:-  Streamline  Panel  Selection. 


Column 

Variable 

Value 

Description 

1-5 

NLXN 

variable 

Number  of  prescribed  panels  for  stream- 
line calculations 

6-10 

NSP  (1) 

variable 

Panel  number  for  streamline  number  1 

11-15 

NSP (2) 

II 

Panel  number  for  streamline  number  2 

16-20 

NSP  (3) . . 

• etc. 

Note: 

Up  to  25  streamlines  may  be  traced  if  only  potential  flow  is 
desired.  Streamline  limit  in  the  boundary  layer  calculation  is  8. 

Note: 

Repeat  Card  34  once  for 
cous  flow  called  for  by 

every  iteration  between  potential  and  vis- 
TRMAX  on  Card  2. 

Card  Set  35;-  This  Card  Set  provides  the  patch  joint  overlap  information 
required  by  the  streamline  routines. 


Card  35A: 

- Number  of 

Joints  to  be 

Input. 

Column 

Variable 

Value 

Description 

1-5 

NJT 

arbitrary 

integer 

Number  of  patch  joints  following 
Set  35B 

Card  Set 

35B:-  Patch 

Joint  List. 

Column 

Variable 

Value 

Description 

1-5 

JOINT (1,1) 

arbitrary 

integer 

Overlapped  patch  number 

6-10 

JOINT (I, 2) 

It 

Side  number  on  overlapped  patch 

11-15 

JOINT(I, 3) 

M 

Overlapping  patch  number 

16-20 

JOINT (I, 4) 

n 

Side  Number  on  overlapping  patch 

Repeat  Card  35B  NJT  times  (I  - 1,  NJT) 
Card  36:-  End  of  information.  6/7/8/9 


79 


2652-79 


